Emergence of ripples on the surface of icicles 
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Natural icicles often exhibit ripples about their circumference which are due to a morphological 
instabihty. We present an experimental study that exposes the mechanism of the instability, using 
laboratory-grown icicles. We find that, contrary to theoretical expectations, icicles grown from pure 
water do not exhibit ripples. The addition of non-ionic surfactants that reduce the surface tension 
does not produce ripples. Instead, ripples emerge on icicles grown from water with dissolved ionic 
impurities. We find that even minute levels of impurity are sufficient to trigger ripples and that the 
growth rate of the ripples increases very weakly with ionic concentration. We outline a proposed 
theory of the instability based on constitutional undercooling. 



Consider the icicle. Its elegant and familiar form is 
the result of a subtle interplay between the solidifica- 
tion dynamics of ice [TJ [2] and the gravity-driven flow 
of the liquid water film flowing over its evolving sur- 
face [3lfT0]. The latent heat released by freezing, which 
controls ice growth, is advected by the water film and 
ultimately carried away by the surrounding sub-zero air, 
which is also flowing [UITT]. Ideal, platonic icicles have 
been predicted to adopt a universal self-similar shape, 
independent of growing conditions [4]. In reality, nat- 
ural [12l HI] and even laboratory-grown [T4lfl6] icicles 
often deviate in interesting ways from the platonic ideal. 
The most striking of these deviations is the formation 
of trains of uniform ripples around the circumference of 
the icicle [S HIOl [TMlG] . Like the overall icicle shape, the 
wavelength of the ripples is remarkably independent of 
the growing conditions. Ripples have been observed to 
move slowly in the direction opposite to the water flow 
during growth [TOl HI [S] . 

Here, we present an experimental study of the physical 
parameters that influence the formation of icicle ripples 
that reveals the essential but hitherto unsuspected role 
of dissolved impurities in their formation. We find that 
ripple growth rates increase roughly logarithmically with 
the concentration of ionic impurities in the feed water. 
Pure water icicles grow no ripples. Dissolved air does 
not, by itself, cause ripples to grow. Reducing the air- 
water surface tension with a strong non-ionic surfactant 
does not produce ripples unless a dissolved ionic impurity 
is also present. These results are in direct contradiction 
to existing theories [5l410j and may find their explanation 
in a theory which incorporates the advection-diffusion of 
both heat and salt, with constitutional undercooling. We 
outline this proposed theory here, which has some simi- 
larity to a recent theory of crenulations on stalactites [T7j . 

In addition to ambient temperature and feed water fiow 
rate [14] . it has previously been shown that water purity 
and air motion are important to the overall icicle mor- 
phology [T^. Icicles grown in still air tend to develop 
multiple tips, while those grown from distilled water are 
significantly smoother and closer to the self-similar shape 
predicted by theory [5] than those grown from untreated 
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FIG. 1: (Color) Images of three icicles grown under identical 
conditions of ambient temperature, water supply rate, and 
nozzle temperature, (a) was made with distilled water only; 
(b) was made with distilled water plus NaCl with a mass 
fraction of (8.0±0.2) x 10"^ (c) was made with distilled water 
plus NaCl with a mass fraction of (1.278 ± 0.002) x 10"^. See 
Supplemental Material at [URL will be inserted by publisher] 
for time lapse videos of the growth process. 

tap water. Tap water icicles show prominent ripples and 
other non-ideal morphologies [16 . Here, we focus on 
comparing icicles grown under the same conditions of 
temperature, water supply rate, and air motion, vary- 
ing only the composition of the feed water. We add small 
quantities of salt to distilled water to achieve levels of im- 
purities similar to tap water; we also investigate both the 
role of dissolved gas and the effect of surfactants which 
reduce the air-water surface tension. We find that the 
presence of small amounts of ionic impurities is crucial 
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FIG. 2: (Color) Spacetime plots showing the time evolution of the ripple topography, viewed from one fixed rotational position, 
for the three icicles shown in Fig. [l] The vertical coordinate is the position on the icicle measured from the root down; the 
horizontal coordinate is time. The edge of the white region shows the downward progress of the tip of the growing icicle. 



to the formation of ripples. 

The apparatus was previously described in Ref. [IS]. 
Icicles were grown below a sharpened wooden support 
suspended inside an insulated, refrigerated box. The sup- 
port was rotated, at a speed of 4 minutes per revolution, 
to encourage axisymmetry and to allow all sides of the 
icicle to be imaged. The air was gently stirred by internal 
fans in the corners of the box. The wall temperature was 
controlled by a commercial bath. The feed water supply 
was delivered by a peristaltic pump to a temperature- 
controlled nozzle that was slightly off-axis from the ro- 
tating support, in order to distribute the water evenly. 
The rotational position of the support was indexed so 
that eight reproducible views of the icicle could be im- 
aged on each rotation. 

We report data from 16 icicles using feed water with 
varying amounts of sodium chloride and surfactant, 
grown under otherwise identical conditions: an ambient 
wall temperature of — 12.2±0.3 °C, a feed water temper- 
ature of 3.0±0.3 °C, and a water supply rate of 2.0±0.1 
g/min. Even from direct inspection, there are obvious 
qualitative differences in the form of the resulting icicles; 
three examples are shown in Fig. [l] Icicles made with 
distilled water have no significant ripples and a smooth, 
carrot-like form that is close to the predicted self-similar 
shape [11[TB]. As the salt concentration is increased, rip- 
ples appear at concentrations similar to, or even below, 
the ionic impurity levels of Toronto tap water. In addi- 
tion to the appearance of the ripples, the overall form 
of the icicle also deviates further from self-similarity and 
axisymmetry with increasing salt concentration |16j . 

We used edge detection to map the topography of the 
icicle in high-resolution digital images, taken at a rate 
of eight images per four-minute icicle rotation. To ex- 
tract the amplitude and wavelength of the ripples as a 
function of time, this data was de-trended using Fourier 
techniques to remove the overall background shape. The 



data acquisition and analysis are discussed further in the 
Supplemental Material at [ URL will be inserted by pub- 
lisher] . 

The resulting de-trended profiles of the three example 
icicles from Fig. [l] are shown in the spacetime plots in 
Fig. [2j Here, the right hand edge of each icicle is shown 
at a particular rotational position. The ripples appear in 
a patchy way and generally grow to a maximum ampli- 
tude early in the lifetime of the icicle. The small residual 
topographic features on the distilled water icicle that sur- 
vive the de-trending filter are consistent with the noise 
floor of the amplitude measurement, while the largest 
features on the saltiest icicle grow to a few mm in am- 
plitude. A trend toward very slow upward motion of 
the ripples is just barely visible in Fig. 2. The overall 
velocity of this motion is consistent with previous obser- 
vations of ripples growing on vertical cylinders [16] and 
inclined troughs |10j . The generally upward crest mo- 
tion is sometimes locally interrupted by the appearance 
of new ripples, ripple mergers, and other dynamics. 

Fig. [3] shows ripple amplitudes vs. time for icicles 
made from various solutions. These were obtained by 
fitting the peak in the spatial Fourier spectrum of the 
de-trended icicle topography, then averaging over one ro- 
tation and over repeated runs under the same conditions. 
We considered ripples only in the uppermost 10 cm of the 
icicle. Ripples farther down the icicle have had less time 
to grow and may grow under saltier conditions due to the 
exclusion of salt by ice formation higher up. In Fig. |3| 
time zero corresponds to the time when the icicle reaches 
10 cm in length. The statistical error bars, which become 
larger for saltier icicles, mainly reflect the progressive loss 
of axisymmetry as the salt concentration increases. Dis- 
tilled water, aerated distilled water, and distilled water 
with Triton X-100 surfactant all show no measurable am- 
plitude increase. The ripple amplitudes for salty icicles 
increase approximately linearly for the flrst 1.5 hours of 
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FIG. 3: (Color) Time-series of the amplitude of ripples on 
icicles made with a variety of salty feed waters. The mass 
fraction of NaCl is shown in the legend. Also shown are pure 
distilled water, aerated distilled water, and distilled water 
with surfactant. For icicles grown from salt water, both the 
ripple growth rate and the maximum amplitude reached in- 
crease with salt concentration. Distilled water without added 
salt shows no significant ripple amplitude growth. The error 
bars indicate the standard error of the mean of rotationally 
averaged data. 

icicle growth. The long-term nonUnear behavior of the 
ripple amplitude is more complex. In many cases, the 
amplitude simply saturates at its maximum value, while 
in others it decays with time. In some cases, ripples con- 
tinue growing slowly until the icicle reaches the maximum 
length allowed by the apparatus. 

Fig.|4]shows some characteristics of the icicle ripples vs. 
the salt concentration of the feed water. In Fig. [4|^a), the 
growth rate of the ripple amplitude is given by the slope 
of Fig. [Sj^a) for the first 1.5 hours. In Fig. |4|^b), we plot 
the error-weighted mean of the ripple wavelength time- 
series after 1.5 hours of growth. As the salt concentration 
is increased, the ripples grow faster. The growth rate in- 
creases only approximately logarithmically with the salt 
concentration. Measurable ripples appear at the remark- 
ably low concentration of 10^^ weight % salt, i.e. only 
10 mg of salt per liter of water. The ripple wavelength 
was 1.04±0.04 cm, independent of the salt concentration. 

We also grew icicles from distilled water that had had 
air bubbled through it for 12 hours, which may therefore 
be assumed to be saturated, or possibly supersaturated, 
with dissolved gases. In fact, all of the solutions used in 
this experiment had long been exposed to air and pre- 
sumably contained some concentration of dissolved gas. 
As shown in Fig. [3]Ja) , ripple formation from deliberately 
aerated distilled water was not significantly different from 
that of untreated distilled water — no significant ripples 



formed. We conclude that dissolved gases alone are in- 
sufficient to trigger the rippling instability. Moreover, 
since the air-water surface tension of the thin water film 
flowing over the icicle surface has featured prominently 
in theories of ripple formation [SHTOl [T7] , we grew icicles 
from distilled water with an added non-ionic surfactant, 
Triton X-100. The results are shown in Fig. [3j Icicles 
grown from distilled water with sufficient surfactant to 
reduce the surface tension by 47% showed no ripples. 
Details on the composition and preparation of the water 
samples is given in the Supplemental Material at [URL 
will be inserted by publisher] . 

Models of icicles 3 -10 assume that their surface is cov- 
ered by a thin film of flowing water which has a parabolic 
shear flow profile with a surface speed U. Film thick- 
nesses h are typically 100 /im or less, with a Reynolds 
number Re = Uh/v ^ 1, where v is the kinematic vis- 
cosity. Any topography present on the ice surface is re- 
fiected in the shape of the water-air surface in a manner 
controlled by the surface tension 7, which appears in the 
free surface boundary conditions in the form of the We- 
ber number We = pU"^ /j ~ 10~^, where p is the density. 
Because the thermal conductivity of air is much less than 
that of water, the temperature difference across the thin 
film AT is much smaller than the total temperature dif- 
ference driving the growth [51 1111 . In our experiments, 
the latter was 12.2° C, while AT < 0.1° C. 

The latent heat released at the ice-water interface is 
diffused and advected through the fiowing water film and 
carried away by the surrounding stirred air. The rele- 
vant dimensionless parameter in the thermal advection- 
diffusion equation for the water film is the Peclet number 
Pch = Uh/n ^ 5, where k is the thermal diffusivity. It 
has been shown that, in addition to diffusion and advec- 
tion, both thermal radiation and water evaporation at 
the water-air interface make significant contributions to 
the heat transport through the surrounding air |11| . 

The rate of ice growth V is contained in the Stefan 
condition pQ [5], which is a boundary condition on the 
temperature gradient VT at the ice-water interface. The 
Stefan boundary condition depends on the Stefan num- 
ber St = L/CpAT 1300, where L is the latent heat 
of fusion and Cp is the heat capacity. St » 1, indicat- 
ing that ice growth is in the quasi-steady limit [T] where 
kinetic undercooling effects should be very small. 

Existing theories of icicle ripples [5-10] have not con- 
sidered the effect of impurities, and included only thermal 
and surface tension effects. They predict ripples on pure 
water icicles, which we do not observe. The linear stabil- 
ity analysis by Ueno [S] found an instability to growing 
ripples that travel against the direction of flow. As is 
usual for such problems, the devil is in the boundary 
conditions. Ueno assumed [10] that the ice-water inter- 
facial temperature obeyed the Stefan condition on VT, 
but that T itself could deviate by 5T away from the equi- 
librium melting temperature T™, which he attributed to 
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FIG. 4: Characteristics of icicle ripples: (a) the amplitude 
growth rate during the first hour; (b) the mean wavelength 
after one hour. Each is plotted against the mass fraction of 
salt in the feed water on a logarithmic scale. Data to the left 
of the break in (a) is for distilled water with no added salt. 



kinetic undercooling [2 [31 [TO] . In this theory, the pre- 
dicted wavenumber kc and the amplitude growth rate <Jr 
of the most unstable mode are given by [8] 
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where g is the acceleration due to gravity. Thus, an in- 
crease in surface tension is predicted to result in a weak 
decrease in the ripple growth rate. Over the range of salt 
concentrations in Fig. [4]^a), however, the air- water sur- 
face tension increases by 0.04% 18 , while the observed 
growth rate increases by a factor of 6. Our results contra- 
dict the 7~i/'^ dependence of Eqn. [l] Furthermore, the 
addition of surfactant, which lowers the air-water surface 
tension, produced no measurable ripples and hence no 
effect on their growth rate. Thus, our results are broadly 
inconsistent with existing theories [5HTO]. 

We now present some general ideas on how to account 
for the ripple growth on salty icicles |19j . In addition 
to the temperature field, the presence of salt in the wa- 
ter film adds another advection-diffusion equation for the 
salt concentration C and a second, much larger Peclet 
number Pcs ~ Uh/D ^ 400, where D is the molecular 
diffusivity of salt. Then, a boundary condition similar to 
the Stefan condition links the production of salt excluded 
by the advancing ice to VC at the ice-water interface. 
This causes a build-up of salt near the interface. This, 
in turn, produces constitutional undercooling [2] , so that 



the temperature at the ice-water interface is given by 
T = Tm — mC, where m is the slope of the liquidus line, 
which accounts for equilibrium freezing point depression. 
This theory ignores kinetic undercooling effects, which 
might still be present for nonzero V. Thus, mC plays 
the role of ST in Ueno's theory [TO], but the physical 
mechanism of the temperature deviation is different, and 
the salt concentration must be non-zero for ripple forma- 
tion. In addition, the detailed effect of the salt depends 
on PCg. 

A natural scaling of the fully double-diffusive problem 
with constitutional undercooling introduces another di- 
mensionless number, M = toCq/AT, where Cq is the con- 
centration of salt in the feed water. M links the bound- 
ary conditions on T and C at the ice-water interface. 
The surprising sensitivity of the ripple instability to very 
small concentrations of salt may be because AT is small, 
only a tiny Cq is enough to make M ^ 1. The absence of 
ripples on pure water icicles presumably reflects the fact 
that their growth rate is small, or possibly even negative, 
for M = 0. It remains to make this idea quantitative. 

There is an imperfect analogy between icicles and sta- 
lactites [TT] [5D]. Stalactites are formed by the deposi- 
tion of CaCOs under similar thin film flow conditions. A 
theory of ripples, known as crenulations, on stalactites 
has recently been proposed |17j . Crenulations also in- 
volve several advecting and diffusing species, with [CO2] 
playing a similar role to the temperature in icicles. The 
dynamics, however, is made more complicated by bulk 
chemical reactions that have no analog in salty icicles. 

In conclusion, we have shown that small amounts of 
ionic impurities are required for the formation of rip- 
ples on icicles and that ripples are not present on ici- 
cles grown from pure water, even with dissolved gases. 
Ripples are not produced by simply reducing the surface 
tension of pure water with surfactants. Ripple growth 
rates show a very weak, roughly logarithmic dependence 
on salt concentration. We have outlined a proposed new 
theory of the rippling instability based on constitutional 
undercooling. 

We thank J. Wettlaufer, J. Neufeld, M. G. Worster and 
R. E. Goldstein for insightful discussions, and C. Ward 
for his assistance with the surface tension measurement. 
This work was supported by the Natural Sciences and 
Engineering Research Council of Canada. 
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Supplementary Information for Physical Review Letters 
Emergence of ripples on the surface of icicles 

Time lapse video of growing rotating icicles 

A supplementary video available on (YouTube) or 
p^'lickr|. This supplementary video shows the evolution 
of the three icicles shown in Fig. 1 in time lapse. Refer to 
the caption of Fig. 1 for the detailed growing conditions. 

Water Samples 

Table |l] shows the measured physical parameters of the 
various water samples used. Toronto tap water is shown 
for comparison. The distilled water was supplied in bulk 
by Canadian Springs . Note that, while it is much 
purer than tap water, it is not near the limits of what pu- 
rity can be achieved. Nevertheless, this level of purity was 
sufficient to suppress the rippling instability. The added 
NaCl was ACS reagent grade, supplied by Sigma- Aldrich. 
The added non-ionic surfactant was Triton X-100, which 
is t-Oct-C6H4-(OCH2CH2)xOH, x = 9-10. The Triton 
X-100 was also supplied by Sigma-Aldrich. 

The conductivity and composition of the samples 
were analyzed by the ANALEST lab, Department of 
Chemistry, University of Toronto. Metal ion concen- 
trations were determined by inductively coupled plasma 
atomic emission spectrometry, using Perkin Elmer Op- 
tima 7300DV. Anions were identified by ion chromatogra- 
phy, using Perkin Elmer Series 200, Alltech Eris lOOOHP, 
AUtech 550 Conductivity Detector, and Phenomenex 
Star Ion 300 Column. Conductivities were measured us- 
ing an Alltech 550 Conductivity Detector. The surface 
tensions of distilled water, tap water, and surfactant so- 
lutions were measured using the capillary tube method at 
room temperature. The surface tensions of salt solutions 
were calculated using Ref. [18]. The results are shown in 
Table D 



Temperature, humidity, and air flow measurements 

The apparatus is a slightly improved version of the 
one discussed in Ref. [16] . The temperature of the inlet 
nozzle was measured with a thermocouple and feedback 
controlled by means of a computer-controlled heater. It 
was maintained at 3.0±0.3 °C to prevent freezing of the 
inlet pipe. The temperature of the feed water had no im- 
portant effect on the icicle growth. The rotating support 
upon which the icicle grew was suspended by stiff wires 
and not actively cooled. All of the heat from the growing 
icicle was transported by advection and diffusion through 
the surrounding air, by radiation, or by evaporation / 



condensation of water [TT]. The latter mechanism de- 
pends on the relative humidity, which was continuously 
monitored during the experiment. For all the data re- 
ported here, the relative humidity was 87±1 %. 

The air inside the icicle-growing box was stirred by 
eight small computer fans: four at the top corners 
pointed straight down, and four at the bottom corners 
pointed straight up. Air motion affects the overall growth 
rates of icicles [14] and is also necessary to suppress the 
tendency of icicles to grow multiple tips |16j . Using a dig- 
ital anemometer (Omega HHF92A), we determined the 
total air flux from each fan to be (1.59±0.05) x 10~^ m'^/s. 
Dividing the volume of the box by eight times this flux 
gives a characteristic time of 0.75 ± 0.05 s — this is an 
estimate of the time it takes for the fans to circulate all 
of the air enclosed in the box once. Thus, the air in the 
box may be regarded as well-stirred, save for a thin vis- 
cous and thermal boundary layer near the icicle and near 
the walls of the box. Thermocouple measurements of the 
air in the well-stirred region show a mean temperature 
of —9.7 ± 0.2°C, a few degrees warmer than the walls 
of the box, which are controlled at — 12.2±0.3 °C by a 
circulating bath. 

Image acquisition and analysis 

To investigate the evolution of icicle ripples, we ana- 
lyzed time series of high-resolution digital images taken 
during the experiment. MATLAB's Sobel algorithm was 
used to detect the left- and righthand edges of the icicle 
in each image. To extract the ripple wavelength and am- 
plitude, the following algorithm was used to de-trend this 
topographic data: the data was first Fourier transformed, 
and the region in the power spectrum corresponding to 
a wavelength between 0.5 - 1.5 cm was set to zero; upon 
back transforming, the resulting ripple-free filtered to- 
pography was smoothed with a moving boxcar average 
with a width of 0.5 cm, giving the overall background 
shape; finally, the background shape was subtracted from 
the original raw edge to obtain the de-trended ripples. 
While the entire de-trended icicle profile is shown in 
Fig. 2, only the topmost 10 cm of the icicle was included 
in the ripple amplitude and wavelength measurements. 

The resolution of the digital images was approximately 
0.018 cm per pixel, so each of the two edges of a 10 cm 
segment of the icicle contributed about 550 edge detec- 
tion measurements. The Fourier power spectrum of the 
de-trended data was fit to a Gaussian peak to extract 
the amplitude and wavelength of the ripples. The data 
was averaged over both edges and one rotation consist- 
ing of 8 unique views, so each data point in Fig. 3 of 
the paper was based on averaging over 8800 edge detec- 
tion measurements. In some cases, the data was further 
averaged over multiple runs under the same conditions. 
The noise floor, which is reached by the residual topog- 
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Solution 


Main impurities 


Concentration 


Conductivity 


Surface tension 


'I'll 1 1 / TXT t T\ 

Distilled water (DW) 


Ca^ 


0.037 mg/L 


2 /iS/cm 


0.072 N/m 






0.026 mg/L 








IN a 


U.Ui4 mg/L 








Ba 


0.003 mg/L 






DW + iriton A-lOO 


iriton X-100 


200 mg/L 


1.7 pS/cm 


0.039 N/m 


DW + NaCl 


NaCl 


80.0 mg/L 


211 ^S/cm 


0.072 N/m 


Toronto tap water 


^ -t-2 

Ca+ 


37.5 mg/L 


419 ^S/cm 


0.071 N/m 




TVT 4- 

Na+ 


10.7 mg/L 










8.82 mg/L 










L66 mg/L 








0-+2 

Sr 


1.22 mg/L 








4-9 


0.46 mg/L 










31.60 mg/L 








cr 


25.11 mg/L 








NO3- 


0.96 mg/L 








p- 


0.41 mg/L 








Zn+^ V+^ Ni+^ 


< 0.2 mg/L 








Ba+^ Fe+^ 









TABLE L Measured compositions, conductivities, and surface tensions of the water samples. Toronto tap water, which produces 
icicles exhibiting prominent ripples, is shown for comparison. 



Physical parameters 






Typical icicle radius 




R 


0.012 m 


Volumetric flow rate 




Q 


3.3 X 10"* m^s 


Typical radial growth speed 


V -- 


= dR/dt 


1.7 X lO"'^ m/s 


Water film thickness 




6.2 X 10"^ m 


Surface speed of water film 


U = 


= gh^ /2i' 


0.011 m/s 


Driving temperature difference 


T 


— Twall 


12.2 °C 


Water film temperature difference 


AT = 


PicoLVh/A 


0.06 °C 


Liquidus slope for salt water 




m 


3.706 °C kg/mol 


Molality of salt 




Co 


0.0219 mol/kg 


Constitutional undercooling 




-mCo 


-0.08 °C 


Latent heat of fusion 




L 


3.34 X 10^ J/kg 


Kinematic viscosity 




V 


1.8 X 10"® mVs 


Density of water 




P 


1.00 X 10^ kg/m^ 


Specific heat of water 




Cp 


4.198 X 10^ J/kgK 
1.36 X 10"^ mVs 


Thermal diffusivity of water 




K 


Molecular diffusivity of salt 




D 


1.6 X 10"'-* mVs 


Surface tension of pure water 




1 


0.0757 J/m^ 


Dimensionless groups 






Ratio of growth to flow speeds 


V 


= V/U 


1.6 X 10"" 


Reynolds number of water fllm 


Re 


= Uh/i' 
= pU^hh 


0.36 


Weber number 


We -- 


9.1 X 10"^ 


Peclet number for heat 


Peh 


= Uh/n 


4.8 


Peclet number for salt 


Pes 


= Uh/D 


410 


Stefan number 


St = 


L/CpAT 


1300 


Constitutional undercooling number 


M = 


mco/AT 


1.3 



TABLE II: Typical quantities and some dimensionless groups for the case shown in Fig. 1(c). Here, A = KpCp is the thermal 
conductivity of water, the values of v, p, 7, k. A, and Cp are from Ref. [18], the value of D is from Ref. [21], and R is taken 
after 2 hours of growth. 



7 



raphy of a distilled water icicle, was approximately 0.005 
cm, or about 1/3 of a pixel. The averaged, best- fit ripple 
amplitudes, shown in Fig. 3, are smaller than the typi- 
cal amplitudes within the patches of de-trended ripples 
shown in Fig. 2, because the fit included Fourier signals 
from regions with no ripples, which reduced the average. 
The errors, shown in Fig. 3, were estimated from the 
standard error of the averaging sample. The whole algo- 
rithm was tested on simulated images which had known 
ripple amplitudes and wavelengths. 

Physical parameters and dimensionless numbers 

In Table |TT] we list the important physical quantities 
likely to be involved in any theory of saline icicle ripples. 
We also estimate the relevant dimensionless groups dis- 
cussed in the paper for a reference icicle - the case shown 
in Fig. 1(c) of the paper. The growth of this icicle is also 
shown in time lapse in the supplementary movie. 

The water film thickness h can be calculated from the 
measured volumetric flow rate of the feed water Q, which 
is held constant. The film thickness can then be used to 
estimate the flow speed U at the surface of the water 
film, treated as having a parabolic flow profile. Most of 
the dimensionless quantities depend on combinations of 
h and U. It important to note that, since the radius 
R of the icicle grows with time, h and U, as well as all 
the dimensionless quantities that depend on them, evolve 
significantly over the lifetime of the experiment. 

To estimate AT, the temperature difference across the 
water film, we use the Stefan condition and treat the film 
as having a perfect parabolic profile so that advection can 
be neglected. Then the heat fiux is purely conductive 
and equal to AAT/h. This must match the rate of heat 
production by ice growth piccLV. We use the measured 
ice growth speed V to estimate AT. 

Since the free surface flow speed U of the parabolic 
profile is much greater than the typical radial ice growth 
speed V, a quasi-static approximation is justified, i.e. 



the steady-state equations for the fluid flow and the 
advection-diffusion of heat and salt are solved over a 
slowly evolving ice topography, approximated as static. 
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